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By extending the response function approach developed in nonlinear optics, we an- 
alytically derive an expression for the non-Markovianity [Laine, et al, Phys. Rev. 
A 81, 062115 (2010)] in the time evolution of a system in contact with a quantum 
mechanical bath, and find a close connection with the directly observable nonlinear 
optical response. The result indicates that memory in the bath induced fluctua- 
tions rather than in the dissipation causes non-Markovianity. Initial correlations 
between states of the system and the bath are shown to be essential for a correct 
understanding of the non-Markovianity. 
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1. Introduction 

For macroscopic systems, the second law of thermodynamics prescribes ever-increasing 
entropy. In fact, decreases of entropy are permitted on short time scales. When the 
dynamics of small quantum systems on short time scales is studied, t he flow of 
inform ation between the environment and the system can be important |Waldraml 
(1985)]. In a microscopic theory, there are three major effects of the environment 



(or bath) on the system. The first two are dissipation, which removes excess energy 
from the system, and fluctuations, which supply energy. These two effects are re- 
lated through the fluctuation dissipation theorem, which assures that the correct 
finite temperature equilibrium state is reached. The third one, which is less well 
known outside the fields of nonlinear optics and NMR, is the presence of entangle- 
ment between system and environment states. This third effect plays a major role 
if the system bath interaction is strong, or if the characteristic time scale of the 
noise induced by the environment is slow compared to typical system time scales. 
It is the origin of a rephasing signal in photon echo and NMR echo measurements. 

The dynamics of a quantum system in contact with a bath is described theo- 
retically by deriving an equation of motion for the reduced density matrix, which 
includes only the system degrees of freedom. It is often assumed that the char- 
acteristic time scale of motion in the environment is much shorter than anything 
that happens in the system. This approximation is convenient, because it allows 
the derivation of a closed equation of motion for the reduced density matrix. In 
particular, no information about the history of the dynamics can be stored in the 
bath, and the equation of motion is local in time. The requirements that the trace 
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of the density matrix must be preserved, and that the diagonal matrix elements 
must be positive in any basis t hen lead to a mast er equation in the Lindblad form 
Breuer fc Petruccione ( 2002 ): Lidar et 'all ( 2001 )]. The dissipation operators that 



appear in this equation can be specified by modeling the environment and its in- 
teraction with the system. Because of the f ast bath approximati on, commonly used 
in association with secular approximations [Fassioli et al. ( 2010[ )]. the master equa- 
tion approach does not properly include the fluctuations caused by the bath. On the 
other hand, stochastic Liouville equation s, which do treat the fluctuations correctly, 
neglect the dissipation Tanimural ( 2006 )]. 

The restriction to an environment with much faster dynamics than the system 
breaks down for many systems studied in ultrafast nonlinear optics and NMR. In 
addition, it is not valid at low temperature, when a correct quantum mechanical 
description of the bath introduces additional time scales determined by the Matsub- 
ara frequencies. Obviously, the dynamics becomes more c omplex in this sit u ation. 
Altho ugh strong system bath interaction can be included Jang et all ( 2008 ) ; Nazir 
(2009)], the key difference with the Lindblad formalism is the presence of memory. 
An environment which is not infinitely fast (compared to typical time scales in the 
system) can store information about the past. This information can subsequently 
flow back into the system, influencing the dynamics. Such memory effects can be 
included in master equations. These can, however, usually not describe a second 
key effect of a slower environment: the presence of correlated super positions of the 
system and the environment in the initial state Suarez et all ( 1992 )]. Because such 
correlations introduce a second source of memory, they cannot in general be ignored. 

The idea that the flow of information from the environment back to the system 
can be used to quanti fy the extent of mem ory i n a non-Markovian quantum process 
has been developed bv lBreuer et al. ( 20091 ) and lLaine et al. ( 2010l ). In a memory less 
situation, two system states that are initially a certain distance apart, will only get 
closer during the time evolution. Therefore, when states are found that grow farther 
apart, the time evolution can be called non-Markovian. By introducing distance 
measures on the space of system states, these ideas can be made precise, resulting 
in a measure for non-Markovianity that depends only on system degrees of freedom. 

Once such theoretical measures have been introduced, the question arises how 
they can be measured in experiment. In prin ciple, quantum st a te tomography 



yields the complete quantum state of the system [Kuah et all (|2007l) ; iMohseni et al 
( 2010l )]. Once this measurement has been performed, any quantity that is a func- 
tional of the reduced density matrix can be calculated. However, this process is 
rather cumbersome and indirect. More straightforward methods to quantify con- 
cepts such as en t anglement and non-M arkovianity from experiments are welcome 
' Xu et all (|2010l ): ICramer et all (|201ll )]. 



Experiments that can achieve this goal are found in the field of nonlinear optics. 
Observables such as the photon echo and two-dimensional optical spectroscopy 
depend on multiple time inter yals, and are sensitive t o memory effects th at extend 
over several of these intervals Mukamel et al. ( 20091 ): Engel et al. ( 2007 )]. 

In this paper, we show how correlations between the system and the environment 
are of critical importance for the non-Markovianity. In particular, for a simple 
model environment the dynamics is completely Markovian if the initial state does 
not include such correlations, while it becomes non-Markovian if correlations are 
allowed to be present. We describe the generation of initial correlations during 
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a preparation time and discuss the close connection with the nonlinear optical 
response. 

2. Trace distance and non-Markovianity 

For a classical stochastic process, the meaning of " Markovian" is clear [Van Kampenl 



(2007)]: the future depends only on the present state, and not on the past. In 
the case of a Gaussian process, its correlation function must be exponential to 
have Markovianit y. In the quantum case, we will define Markovianity following 



Laine et "all (l2010h . although other approaches have been proposed as well (Rivas et al. 



( 2010h ]. As explained in the Introduction, the definition is based on the distance 



between a pair of quantum states. 

A convenient measure is given by the trace distance D between two density 
matrices g A and g B , which is defined as 



D(g A , q b ) = l -Tr s ^{Q A -Q B )\ (2.1) 



2 

or half the sum of the square root of the eigenvalues of (g A — g B ) 2 - The subscript 
S in the trace indicates that it is taken over system degrees of freedom, in contrast 
with trace operations over the environment which we will encounter later. 

In an ergodic system, any initial state will evolve in time until it reaches a 
single well-defined equilibrium. If there is no memory in the bath, the dynamics 
can only bring the system closer to equilibrium. Because no information can flow 
from the environment to the system, the distance between a pair of initial states will 
decrease with time. This is the case in memory-less approaches such as the Lindblad 
master equation. Memory in the bath means that the bath stores information about 
the system at a previous point in time, which affects the dynamics. The extra 
information opens the possibility of temporary time evolution in the unnatural 
direction. This suggest that non-Markovianity can be measured by st udying how 



much two states move away from each other. In the definition given bv lLaine et al 



(2010), this quantity is studied by defining the change in the trace distance 



a(t; /(G), g B (0)) = ±D( Q A (t), g B (t)), (2.2) 
at 

which is integrated over time to define the non-Markovianity 

M(G) = max / dta(t;g A {0),g B (0)). (2.3) 
Jcr>a 

The maximum is taken over all combinations of initial states g A (0) and g B (0). This 
measure is based on the distinguishability of the two trajectories. 

Although the trace distance can be defined for two density matrices of any 
shape, it takes a particularly simple form for a two-level system. In this case, the 
matrix elements of a density matrix are written in a given basis as 

9={ Q1 : / 12 V (2.4) 

The trace distance between g A and g B is found as 



D(g A ,g B ) = J(g A - gf^ + \g A - ef 2 | 2 . (2.5) 
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(a) Initial correlations 

Although the (statistical) state of the system's degrees of freedom at any point in 
time is completely described by the reduced density matrix, this is nevertheless not 
the complete story. Durin g the t ime evolution, the system gets entangled with the 
environment [Lopez et al. (2008)]. If the time scale of the environment is not very 



short, this will influence the state of the system at later points in time. Therefore, 
the presence of classical or quantum mechanical correlations between system and 
environment affects dynamic measures like the non-Markovianity, although they 
are not explicitly present in the reduced description. Although correlations are 
automatically produced during the time evolution, leading to non-Markovianity, 
they can also be present in the initial state. These initial correlations contribute to 
the non-Markovianity as well and should be included in a proper description. To 
see this explicitly, we denote the complete density matrix including all system as 
well as bath degrees of freedom as R(t). Its matrix elements in the system subspace 
are still operators on the bath degrees of freedom. Because the system and bath 
taken together form a normal quantum system, the complete density matrix evolves 
coherently in time, as dictated by the complete Hamiltonian H. We can define a 
propagator G which propagates the density matrix as R(t) — G(t—to)R{to), which is 
given by G(t-t ) = exp(-ii/ x (t -t )/H). The notation H X A = [H,A] denotes the 
commutator. The reduced density matrix, which operates only on the Hilbert space 
of the system, is found by taking the partial trace over the bath, g(t) = ti sR(t)- For 
a factorized initial state, R(to) = g(to)RB(to), the time evolution can be written as 
a dynamical map g(t) = $>(t;to)g(to)- 

However, in the case of a slow environment, it is not clear why initial correlations 
between system and bath states can be neglected and the factorization assumption 
may break down. This means that the complete density matrix cannot be written in 
the form i?(0) = q(0)Rb(0), where i?s(0) is a density matrix in the Hilbert space of 
the bath. Instead, each matrix element in the system space may depend on the bath 
in its own way. The difference between uncorrelated and fully correlated equilibrium 
density ma trices, which are giv en by exp(— /3if)/Trcxp(— (3H), has been studied 
recently by [Sm irne et all hOld ) and is readily observable in the optical response 



Uchivama fc Aiharal (|2010D ]. 

To study the effect of more general initial correlations, we introduce a prepa- 
ration time. This method allows us to interpolate between an uncorrelated state 
and the properly correlated equilibrium. An initially uncorrelated state is allowed 
to evolve for a time t\, during which correlations are formed. The dynamics of the 
thus obtained correlated state at time zero is then followed during a time t%. For 
t\ = 0, correlations between system and bath are absent, while for a long enough 
preparation time maximum correlation is reached. 



3. Model 

The coherent time evolution of the system is given by a Hamiltonian H$ . We employ 
a commonly used model for the environment that includes the complete quantum 
mechanical behaviour of bath modes, yet is flexible enough to be solved to a certain 
degree. In this model, the bath modes are harmonic oscillators, which couple linearly 
to the system. The Hamiltonian for the bath and its coupling to the system is given 
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(3.1) 



Here, a indexes the bath modes, which have coordinates x a , momenta p a and 
masses m a . V denotes any operator on the Hilbcrt space of the system, which cou- 
ples to the bath modes with strength g a . All necessary information about the system 

2 

bath interaction is contained in the spectral density J(lo) = ^ a S(lu — uj a ) 

and the temperature T. The correlation function can be written as the inverse 
Fourier transform of the spectral density as 



1 f°° 

Lit) = — I doj J(uj) (coth/37ia;/2 coswt — isinwt) . 

1" JO 



(3.2) 



Its real part corresponds to the fluctuations, which are a function of the inverse 
temperature /? = 1/fcsT (fcs is the Boltzmann constant), whereas the imaginary 
part is the dissipation. Because the statistics for linear coupling to a harmonic 
bath is the same as for a Gaussian process, multi-point correlations functions are 
redundant (they can be evaluated using Wick's theorem) , and are not necessary for 
the calculation of the propagator. In the classical limit, the bath can be modeled 
by a stochastic process. The time evolution is t hen given by a sto chastic Liouville 



equation, which can include initial correlations Ban et al 



for the calculation of the nonlinear response [Jansen et al 



(20106)] 
(|2004h ]. 



and is suitable 



To simplify the analytical treatment, we will describe the situation where the 
system-bath interaction commutes with the system Hamiltonian, [V, Hs] = 0, 
such that the exact dynamics be c omes second orde r in the sys t em-bat h interac- 
tion [ishizaki fc Tanimural (|2008h ; iNan et all (|2009t ): iBan et all (|2010ah ]. For an 
overdamped Brownian oscillator spectral density, the non-commuting c ase can be 
handl e d efficiently using the hierarchy of equations of motion approach Tanimura 
(|2006l ): IShi efaZI (|2010h ]. 

Although the system Hamiltonian can be chosen freely, we will for definiteness 
focus on a two-level system. In the basis of its eigenstates, the system Hamiltonian 
is diagonal, with matrix elements and e. The system-bath interaction causes 
dephasing in the excited state, and Hsb has matrix elements and Se(X), where 
the fact that Se is an operator on the bath degrees of freedom is indicated explicitly 
by the notation (X). 

In a linear response experiment, the system is brought out of equilibrium by 
an external pulse, and the subsequent time evolution is probed. Non-Markovianity 
during the evolution time can occur in two ways. Firstly, it can be caused by 
memory in the system bath interaction during the evolution time. A second source of 
non-Markovianity are initial correlations between the system and the environment, 
which are present at the time the impulsive force interacts with the system. Such 
correlations can be studied in detail using nonlinear experiments, involving multiple 
pulses. 
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4. Results 

(a) Trace distance as a function of a single time 

When a system is initially in a factorized state, the only source of non-Markovianity 
is the buildup of system bath correlations during the time evolution. Suppose 
that two initial density matrices are given by R A (Q) and R B (Q). We make the 
usual assumption (which we want to relax later) that the system can be separated 
from the bath, and that the bath is in thermal equilibrium. The complete density 
matrix is then written as the direct product of a system part and a bath part, 
R A (0) — £> (0) exp(— (3Hb)/T^b exp(— (3Hb), where the reduced density matrix is 
g A {0) — Trsi? A (0). Similar relations are written for R B (0). Such factorized ini- 
tial conditions are typically found in electronic resonant spectroscopy, where the 
thermal energy is much smaller than the electronic excitation energy. The equilib- 
rium density matrix then only contains population in the ground state, given by 
Qcq = From this state, one can create any factorized initial state by applying 

an impulsive external interaction g A / B (fS) = Ug cq , where U denotes a Liouville 
operator. 

The time evolution of the complete density matrix is given by coherent evolution 

R A (t) = exp(-iH x t/h)R A (Q) = exp(-i( J ff| + H* + Hg B )t/h)R A (0). (4.1) 

We assume that the system Hamiltonian commutes with the system-bath inter- 
action. In the interaction picture with respect to the bath Hamiltonian, the time 
evolution of the reduced density matrix then becomes 

g A (t)= C M-^t/h)( C x P+ (-^l\rH^ B (r)))g A (0), (4.2) 

where (• • • ) = Tr^ • • • exp(— /3Hb)/Ttb exp(— (3Hb) and exp + denotes the time 
ordered exponential. The evaluation of the average over the bath is a standard 
calculation, giving 

(exp + (-^ t drff s x i(r))) = J* j* dt"{H* B (t')H* B (t")}), (4.3) 

which for the coherences reduces to exp(— g(t)), while the populations are constant 
in time. The dephasing function is given by 

9(t) = ± I fdt'J t dt"L(t"), (4.4) 

with the correlation functi on L(t) giv e n by F ourier transforms of the spectral density 
according to equation 13.21 [Mukamell (|l995l) ]. 

For a two-level system, the trace distance between two density matrices which 
are initially prepared as g (0) and g B (0) can be readily evaluated using equation 
12.51 It is found to be 

D(g A (t),g B (t)) = yViW - ^(O)) 2 + K 2 (0) - e ? 2 (0)\ 2 exp(-2Re.g(i)). (4.5) 

Because we study the pure dephasing case, the populations g A ^ and g^ are constant 
in time, while the coherences evolve according to the dephasing function g(t). While 
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the dephasing function contains an imaginary (dissipative) part, which causes a 
time-dependent shift in the effective frequency, only the real part appears in the 
trace distance. If two state are prepared at time zero with initial populations, the 
trace distance simplifies to 

D(g A (t),g B (t)) = D(g A (0) )Q B (0))exp(-Reg(t)). (4.6) 

In this case, the trace distance is directly related to the dephasing function. It 
is important to notice that the non-Markovianity only depends on the fluctuation 
part of the bath contribution, represented by the real part of g(t). Thus, one cannot 
reveal this effect from Lindblad-like quantum master equations, which only includes 
the dissipative part of bath contribution properly. However, stochastic Liouville 
equations may be useful for the study of non-Markovianity. 

We are now in a position to analyse the conditions for which the dynamics 
is non-Markovian. According to the definition of the non-Markovianity, the dy- 
namics is non-Markovian only if the trace distance between two density matri- 
ces increases with time. The time derivative of the trace distance is found to be 
D(t) = —Reg(t) exp(— Re g(t)). Because the exponential of a real number is always 
positive, the time derivative can be positive only if Re<7(i) < 0. From the defini- 
tion in equation 14.41 the time derivative is g(t) — J Q dri(r)/?i 2 . We see that the 
trace distance can only increase if the real part of the correlation function is neg- 
ative, and sufficiently negative. Although the relation between the trace distance 
and the dephasing function is more complex in the general case where g A x ^ gf 1: 
the populations do not influence the question whether the dynamics is Markovian. 
Non-Markovian time evolution is found if the trace distance increases at a certain 
point in time. Because the trace distance is a positive quantity, its derivative is given 
by a positive constant times — Reg(t) also in the case of different populations. The 
previous analysis therefore applies, even though the value of the non-Markovianity 
will be different. 

However, because we started from initial states where the system and bath are 
factorized, this treatment does not include initial correlations between the system 
and the bath. To study their effect, we next study the non-Markovianity after an 
initial preparation time. 

(b) Trace distance as a function of two times 

To include initial correlations, we consider a preparation time. Starting from 
an state that factorizes into system and bath parts, which can be created as g° — 
Ug cq — U\l)(l \ in optical experiments, the sample evolves during a time t%. During 
this time correlations between the system and the bath form. An impulsive external 
force U' is then applied to the system, after which time evolution takes place during 
an interval £2 ■ The time variables are illustrated in figure [TJ The non-Markovianity 
during the time £2 can now be caused by two effects: correlations that build up 
during £2, as well as initial correlations present at the moment the external force 
interacts with the system, which are the result of the preparation. 

The density matrix after evolution during two times is given by g(ti,t2) — 
TrBG(t2)U'G(—ti)g°, where we assume factorized conditions at time —t\. U' de- 
notes a Liouville operator that models the second impulsive external force. The 
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matrix product in the system Liouville space can be worked out explicitly by choos- 
ing a basis. We order the basis states as |1)(1|, |1)(2|, |2)(1|, |2)(2| and denote the 

Assuming that [Hg, Hsb] = 0, we find 



matrix elements of U' in this basis as T]L kl 



g(h,t 2 ) = Tr B 



/ #11,11 

G #12,11 
C2 #21,11 
\ #22,11 



C*TV 

S.1 u ll,12 
CICl #12,12 
CaCl #21,12 
Cl #22,12 



& #11,21 

CI Cl #12,21 
C2C1 #21,21 
Cl #21,21 



TP 

u ll,22 
Cl #12,22 
C2 #21,22 
IT' 

u 21.22 



(4.7) 



where 



Ci 

C2 



exp(-ieii/7i)exp + (-- J dr5e(X(r))), 
i f t2 

exp(-iet 2 /h) exp,(-- / dT<5e(X(r))) 
^ Jo 



(4.8) 



are still operators on the bath degrees of freedom. 

The trace over the bath degrees of freedom can now be calculated analytically, 
using cumulant expansion or path integral methods. It results in dephasing func- 
tions, which depend only on a single time when the average over either £1 or £2 is 
taken, but explicitly on both times for the average of products of two £ functions. 
Because we are interested in the effect of initial correlations, these terms, which 
cannot be factorized into separate contributions depending on t\ and t 2 only, are 
the most relevant to our treatment. They contain the effect of memory that extends 
over the externally applied force. From equation 14.71 we see that these interesting 
terms multiply the matrix elements of the external force that operate on the co- 
herences £>i2 and 021- There are four such terms, two which leave the coherence 
unchanged, and two which interchange the two coherences. To focus clearly on the 
effect of initial correlations, we choose an operation that flips the coherence, while 
leaving the populations unaffected. Such a force is given by a Liouville operator 



with matrix elements U'n n = U 22 



U( 



1, and all other elements 



zero. Writing out the matrix elements in equation 14.71 the density matrix is then 
given by 

g(h,t 2 ) = (eSi,SCi^i,CaCrrf 2 .&)- ( 4 -9) 

Starting from two density matrices g (0) and g B (0), with equal initial popu- 
lations, the trace distance between them evolves in time as D(g A (t 2 ), g B {t 2 )) — 
D(g A (—ti), g B (~ti))\Qd\. As in the previous case of a single time interval, the re- 
striction to equal populations only changes the value of the trace distance, but not 
the question whether the dynamics is Markovian. Using the cumulant expansion, 
we find 

IC2C1I = exp [-2Re<?(ii) - 2Reg{t 2 ) + Reg{h + t 2 )] . (4.10) 

This expression enables the straightforward evaluation of the trace distance and the 
non-Markovianity for any spectral density. The term g(t\ + t 2 ) indicates the effect 
of initial correlations present at the time of interaction with the impulsive force 
U'. Such correlations, which extend across the the excitation, cannot be treated 
by the conventional reduced equation of motion approach, which includes Redfield 
and Lindblad equatio ns. This has been poin t ed ou t in the calculations of nonlinear 
optical observables bv llshizaki fc Tanimura ( 2008 '). 
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(c) Two-level system with overdamped bath 



As a simple example, which allows a more detailed analytical treatment, we 
will discuss the case of an overdamped Brownian oscillator. The spect ral density 



is given by J( o j) = 2 nuj^/(oj 2 + 7 2 ), which gives the dephasing function Tanimura 
(ll990h : IWeisJ d2008h ] 



g(t) = -icot(^87/2)(exp(-7*)+7t-l) + 

7 

- i(7j/7)(exp(-7t)+7i-l). 



4?77 exp(— v n t) + v n t - 1 



«r/7 \ - 

h/3 2-~> 



v n{vl -7 2 ) 



(4.11) 



For this model, the imaginary (dissipative) part of the correlation and dephasing 
functions depends only on the single time scale 7. As can be seen from the first 
line of equation 14.111 the real (fluctuation) part includes additional time scales 
dictated by the Matsubara frequencies v n = 2im/ fih. While the imaginary part of 
g(t) becomes constant due to the Ohmic nature of J(ui) for 7 — > 00, the real part 
is time dependent as long as j3 is small. This indicates that the fluctuation part of 
the bath noise cannot be delta correlated, even if this approximation is valid for 
the dissipation part. At high temperature compared to the time scale of the bath, 
Hf3"f/2 <C 1, these quantum fluctuation terms can be dropped, and the dephasing 
function simplifies 

g(t) = (2n/(3h 1 2 )(cM-lt) +lt-l)- (exp(- 7 £) +yt- 1). (4.12) 

In the high temperature case given above, the time derivative of the dephasing func- 
tion is clearly always positive, and, consequently, the non-Markovianity vanishes for 
a single time interval. 

However, if we account for initial correlations by allowing them to form during 
a preparation time, the dynamics can become non-Markovian. This can be seen ex- 
plicitly using the overdamped Brownian oscillator spectral density. For t\ = 0, the 
trace distance varies with time as exp(— Re 5 (£2))) and we recover the result found 
earlier. Because Reg(t) > for all times, the trace distance is strictly decreas- 
ing. Thus, the dynamics is Markovian for an exponential correlation function, in 
agreement with the classical definition. On the other hand, if we allow system bath 
correlations to form during the preparation time £1, the trace distance can increase, 
and the measure for non-Markovianity is nonzero. This effect is shown in figure [2j 
where we compare a factorized initial state, corresponding to t\ = 0, to a state 
that contains correlations, created by setting t\ ^ 0. For a nonzero preparation 
time, the trace distance increases during a certain time interval, which shows that 
the dynamics is non-Markovian. It is clear that memory in the bath that extends 
over the pulse is crucial for this effect. The time evolution of the system after the 
application of the pulse is influenced by its state before the pulse, as can be seen 
from the presence of the g(t\ + £2) term. For a bath that contains memory, this 
term does not factorize into functions of £1 and £2 only. The memory effect can be 
seen clearly by looking at the trace distance as a function of both times, as plotted 
in figure [3] 
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(d) Nonlinear optical response functions 

The time evolution of a quantum system during two intervals, separated by 
an external impulsive force, is closely related to nonlinear optical experiments. In 
these experiments, an initial pulse excites the system out of the ground state. Cor- 
relations between the system and the bath form during the following propagation 
time. After applying another pulse, the effect of these correlations can be observed. 
As mentioned before, the initial state g° can be prepared by applying a pulse on 
the equilibrium state, g° = U g eq . In the case of optical experiments, the excitation 
energy is typically much larger than the thermal energy, and the equilibrium dis- 
tribution only contains population in the ground state. The optical held couples to 
the dipole of the system, described by the dipole operator p, = fl(\l)(2\ + |2)(1|). If 
we choose the operator U as the commutator of the dipole with the density matrix, 
U g = [ft, g], and furthermore set U' — U 2 , the time evolution operators evaluated in 
the previous section correspond t o the standard thi r d-order nonlinear response fun c- 
tions with zero population time Mukamel (|l995lh Tanimura fc Okumura (1996)]. 



The observable in nonlinear optical experiments is the trace of the dipole operator 
multiplied with the density matrix, Trgfigftx, t%) = TrsP^BG(t2)U'G(—ti)Ug oq . 
By choosing the wave vectors of the incident pulses properly, it is possible to select 
pathways that are sensitive directly to the coherence flip described in t he pr evious 



section. The resulting photon echo signal is given by [Nibbering etali (|l99lh 



R(t) = exp [-2ff(ti) - 2g(t 2 ) + g(t x + t 2 )} . (4.13) 

Using heterodyne detection, both the real and imaginary parts of this response func- 
tion are observable, while homodyne detection directly yields the absolute value. 
The connection with the non-Markovianity in the previous section is immediately 
clear: the photon echo is sensitive to exactly the memory effects that are responsi- 
ble for non-Markovian dynamics. A photon echo experiment can be used to prepare 
a state in which the system and environment are correlated, and to subsequently 
probe the time evolution. Plotting the thus obtained response function directly an- 
swers the question whether th e dyna mics is Markovian or not, according to the 
definition given bv lLaine et 



5. Conclusion 

We have studied the non-Markovianity in quantum mechanical time evolution. This 
concept of Markovianity can be made precise by looking for states that become 
more distinguishable during time evolution. If such states are present, the process is 
clearly non- Markovian, which i s the basic idea of the measure for non-Markovianity 
proposed by Laine et al. ( 2010f ). Non-Markovian time evolution corresponds to the 



presence of memory effects. Only the (temperature dependent) bath induced fluc- 
tuations, and not the dissipation, enter the non-Markovianity. We have treated the 
dynamics without making Markovian or rotating wave approximations, and thereby 
fully included the correlations between system and bath states, which influence the 
dynamics at a later point in time. Not only the correlations that are formed during 
the evolution, but also those present in the initial state can cause memory effects. 
By forming a correlated initial state during a preparation time, this effect can be 
studied for more general initial states than the equilibrium with respect to the com- 
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plete Hamiltonian. We have shown that a process that is Markovian without initial 
correlations can become non-Markovian when such correlations are present. Con- 
ventional master equations, which cannot include the preservation of system bath 
entanglement across a pulse cannot be used to analyse this situation. Clearly, com- 
monly used approximations such as a delta correlated bath or secular system-bath 
interaction don't hold either. Because the procedure of preparing correlations dur- 
ing an initial time, and subsequently measuring their effect following an external 
impulse, the non-Markovianity is directly observable in nonlinear optical experi- 
ments such as the photon echo. Future work should consider the three time photon 
echo, and the closely related two-dimensional optical spectra. In these experiments, 
population dynamics can be studied during a waiting time, allowing for more gen- 
eral measures of non-Markovian time evolution. Generalizations of the current work 
to more general system Hamiltonians, multiple baths, the case where the system 
Hamiltonian and the system-bath interaction do not comm ute and low tempera- 



ture a r e possible using t he hierarchy of equati o ns of motion [Ishizaki fc Tanimura 
(l2005h : IShi eTaH (l2010h : fPiikstra fc Tanimural d2010h ]. 
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Figure 1. Schematic drawing indicating the time variables and the state of the reduced 

density matrix. 
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Figure 2. Trace distance as a function of time for preparation times t\ — (dashed line) 
and ti = 1 (solid line). The environment is modeled as an overdamped Brownian oscillator 
at high temperature with parameters f3tvy = 0.5 and /3hrj = 1.0. The dotted lines show 
the same quantities calculated with 100 low temperature correction terms. Increase of the 
trace distance with time shows non-Markovian evolution. 
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Figure 3. Trace distance as a function of two times, a preparation time t\ and a detection 
time t2- Parameters for the environment are the same as in figure [2] Correlations between 
the time evolution during ti and £2 are clearly present. 
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